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promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


Hydraulic  wicket  gates  are  modem  appurtenances  for  flow  control  at  locks- 
and-dam  systems  while  maintaining  the  navigational  characteristics  of  a  pass. 
Stmcturally,  these  gates  are  a  group  of  edge-hinged  steel  panels  acting  in  parallel 
and  independentiy  driven  by  remote-controlled  hydraulic  actuators.  The  gates 
operate  from  a  horizontal  open  position  to  an  engineered  elevation  in  the  closed- 
gate  condition.  Due  to  complex  surrounding  flow,  individual  gates  are  subjected 
to  hydraulic  forces  which  locally  compound  the  steady-state  component  of  the 
induced  stresses  with  random  fluctuations. 

The  Olmsted  locks-and-dam  system  will  use  220  hydraulic  wickets  for 
navigational  flow  control  of  the  Ohio  River  at  Paducah,  KY  (March  and  Elder 
1992).  This  report  describes  the  experimental  analysis  to  define  a  mathematical 
model  to  estimate  the  hydraulic  forces  on  a  prototypical  Olmsted  wicket  gate 
under  upstream  quasi-laminar  flow.  In  the  mathematical  model,  water  is 
assumed  an  ideal  homogeneous,  incompressible  and  inviscid  fluid,  and  its  flow  is 
taken  as  irrotatirmal.  Quasi-laminar  flow  is  defined  as  flow  associated  with  small 
to  moderate  velocities  (Reynolds’  number  smaller  than  approximately  2000 
(Tritton  1982)). 

The  hydraulic  environment  in  which  water  flows  through  a  three-gate 
ape^re  and  subsequently  impinges  on  a  single  central  gate  in  closed  position 
(60°  eleva-tion)  was  identified  from  experimentation  on  a  corresponding 
1 /25-scale  physical  model  at  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  as  one  of  the  critical  conditions  to  be  considered  in  the  gate 
design. 

Figures  1  through  3  illustrate  the  setup  for  the  experiments.  Figure  1  shows 
the  diy  physical  model.  Figure  2  illustrates  an  instance  of  quasi-laminar  flow 
through  a  three-gate  aperture,  and  Figure  3  illustrates  the  quasi-laminar  hydraulic 
environment  acting  on  a  wicket  gate  at  45  °  elevation. 
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2  Theoretical  Development 


Eccentricity  of  the  Conics 

The  notion  of  eccentricity  (or  flatness)  of  the  conics  is  briefly  reviewed  and 
summarized  as  the  starting  step  of  the  development  (Kreyszig  1988). 

Ellipse 

Figure  4a  shows  a  general  ellipse  with  equation; 


(1) 


where  5  and  t  are  the  major  and  minor  semiaxes,  respectively. 

Points  P  on  the  ellipse  have  the  property  that  the  sum  of  the  distances  from  P 
to  the  foci  (Fj  and  Fj)  is  constant.  The  value  of  this  constant  is  apparent  when  the 
end  A  of  the  major  semiaxis  is  considered; 

dy  +  flfj  ~  2  5  (2) 

Consideration  of  end  B  of  the  minor  semiaxis  leads  to  an  expression  for  the 
focal  distance  c; 


-is 


C  = 


(3) 


where  s>  c. 

The  eccentricity  of  the  ellipse  is  defined  as; 


e  =  —  =  COS0  ^  1 
s 


(4) 


The  completely  nonflat  ellipse  is  a  circle,  for  which  e  =  0,  and  the  completely 
flat  ellipse  is  the  focal  segment  FjFj ,  for  which  e-\. 
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Figure  1.  Physical  model  of  Olmsted  wicket  gate  (Scale  1:25) 
Chapter  2  Theoretical  Development 
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Rgure  4.  Eccentricity  of  conics:  (a)  Ellipse,  (b)  Hyperbola,  and  (c)  Summary 
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Hyperbola 


Figure  4b  shows  a  general  hyperbola  with  equation: 

£!-z!  =  i 


(5) 


where  s  and  t  are  the  transverse  and  the  conjugate  semiaxes,  respectively. 

Points  P  on  the  hyperbola  have  the  property  that  the  absolute  value  of  the 
difference  between  the  distances  from  P  to  the  foci  (F,  and  Fj)  is  constant.  The 
value  of  the  constant  is  apparent  when  the  end  A  of  the  transverse  semiaxis  is 
considered; 


1^1  “  41  =  25 

Consideration  of  the  end  B  of  the  conjugate  semiaxis  leads  to  an  expression 
for  the  focal  distance  c: 


(6) 


{s 


where  Os. 

The  eccentricity  (flatness)  of  the  hyperbola  is  defined  as: 


(7) 


e  =  —  =  sec0  ^  1 

5 


(8) 


The  completely  flat  hyperbola  is  the  discontinuous  real  axis  given  by  |x|  s  c, 
for  which  e  =  1,  and  the  completely  upright  hyperbola  tends  to  the  y-axis,  for 
which  e-oo. 


Figure  4c  summarizes  the  classification  of  the  conics  according  to  their 
eccentricity. 


Laminar  Flow  Through  Fundamental 
Aperture  at  Boundary  of  Half-Plane 

Consider  a  plane  (of  infinite  extension)  divided  in  two  by  a  straight  and  rigid 
boundary  with  an  aperture  of  size  2  length  units  centered  at  the  origin  of  cartesian 
coordinates,  as  shovm  in  Figure  5.  The  two-dimensional  laminar  flow  through 
this  fiindam^tal  aperture  is  conveniently  described  in  terms  of  potential  theory 
(de  Bejar  and  Hall  1994). 
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The  complex  potential  in  the  plane  of  flow  is  defined  as: 

where 

y)  =  potential  function  (Lamb  1945,  Prandtl  and  Tietjens  1957) 
y) ""  corresponding  streamline  function 
Z=  general  complex  variable 
i  =  imaginary  unit 

The  following  complex  potential  proves  to  satisfy  all  boundary  conditions  of 
the  problem  and  therefore  renders  its  solution: 


F(Z)  =  cosh'^Z  =  ln(Z  +  -  1)  (10) 

Putting  Equation  9  into  Equation  10,  a  system  of  two  real  equations  is 
obtained: 

e^*-cos2;|f  -  2e*-(xcos\|;  -ysinr|;)  +  1=0 
e ^ •  sin 2i|r  -  2e‘*’'(Arsini|;  +^cosi|;)  =  0 

Elimination  of  <|)(x,  y)  in  Equations  1 1  leads  to  the  following  geometrical 
locus,  after  some  algebraic  manipulations: 


cos^  ij;  sin^  t|; 


which  defines  a  family  of  confocal  hyperbolas:  the  streamlines.  In  fact,  the  focal 
distance  for  each  individual  hyperbola  is  given  by  Equation  7: 

c  =  Vcos^t  +  sin^ij;  =  1  (independent  o/f)  (13 

Therefore,  the  ends  of  the  aperture  are  the  foci  for  the  whole  family  of  hyper¬ 
bolic  streamlines.  The  eccentricity  of  each  streamline  is  determined  by  the 
specific  value  of  rfr  and  is  given  by  Equation  8  as  e  =  sec  ijr  ^  1. 

The  corresponding  equipotential  lines  may  be  obtained  by  eliminating 
(x,  y)  in  Equations  1 1,  leading  to: 
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(14) 


cosh^  <|>  sinh^ 

which  defines  a  family  of  confocal  ellipses. 

The  focal  distance  for  each  individual  ellipse  is  given  by  Equation  3: 


c  =  )/cosl?^-sinh^  =  1  {independent  of  4>)  (15) 

Therefore,  the  ends  of  the  aperture  are  also  the  foci  for  the  whole  family  of 
elliptical  equipotential  lines  which,  together  with  the  streamlines,  define  a 
network  of  confocal  orthogonal  conics.  The  eccentricity  of  each  equipotential 
line  is  determined  by  the  specific  value  of  <J)  and  is  given  by  Equation  4  as 
e  =  sech<j)  1. 


Laminar  Flow  Through  Finite  Aperture 
at  Boundary  of  Half-Plane 

Similarly,  considering  the  two-dimensional  laminar  flow  through  a  finite 
aperture  of  size  2a  within  a  straight  and  rigid  boundary  dividing  in  two  the  plane 
of  the  flow  leads  to  a  network  of  confocal  orthogonal  elliptic  equipotential  lines 
and  hyperbolic  streamlines,  as  Figure  6  shows. 

The  conforming  complex  potential  is  now  given  by. 

F{Z)  =  Kj- cosh'  *  ^  •  In  (^  ^  ~  °-)  (16) 

where  Kj  is  a  constant  with  appropriate  dimensions  [length  ^/time]  to  be 
determined  from  the  actual  flow  velocity  measured  at  a  reference  point  during  the 
experiments  on  the  physical  model. 

With  the  definitions  <l>'  (x,  y)  =  <j)  (x,  y)/  Kj  and  i|f'  (x,  y)  =  i|f  (x,  y)/  Kj ,  the 
family  of  streamlines  is  given  by  the  hyperbolas: 


a^-cos^ijf'^  a^*sin^i|r^ 


and  the  family  of  equipotential  lines  is  given  by  the  ellipses: 


a^’cosh^^'^  a^*sinh^<|)^ 


=  1 


(18) 


10 


Chapter  2  Theoretical  Development 


11 


Figure  6.  Two-dimensional  laminar  flow  through  finite  aperture 


The  common  foci  for  this  network  of  conics  are  located  at  the  ends  of  the 
finite  aperture  (i.e.,  x  -  ±a).  The  eccentricity  of  the  streamlines  is  given  by 
e  =  sec  i|;'  s  1,  and  the  eccentricity  of  the  equipotential  lines  is  given  by 
e  =  sech<j>'  ^  1. 


Mass  Flux  Through  Finite  Aperture 

The  mass  flux  through  a  finite  aperture  of  size  2a  and  unit  depth  equals,  by 
continuity  arguments,  the  mass  flux  through  the  upper  half  of  an  arbitrarily 
selected  equipotential  line  of  length,  say,  A  (Figure  6).  This  flux  is  given  by: 

^  =  2pj’  |{o)}|  dK  (19) 

A 

where  (w)  =  u  +  i*v  is  the  complex  velocity  field. 

The  complex  conjugate  velocity  field  is  given  by: 

{CO}*  =  F'{Z)  =  [cosh-'^l  (20) 

aZ  a 

Equating  real  and  imaginary  parts  of  Equation  20  yields  a  system  of  two  real 
equations  for  the  components  u  and  v  of  the  velocity  field: 


(21a) 

(21b) 


Since  j  {a>}  |^  =  +  v^,  the  integrand  in  Equation  19  is  obtained  as: 


{(x^  -  -  a^f  +  4x^y^}^ 

and  the  arc  length  may  be  obtained  by  implicit  differentiation  of  Equation  18  as: 


12 
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Substitution  in  Equation  19  gives  an  expression  for  the  flux: 


a  cosh 

q  =  2p -Kj- seek  j 


a^cosh"*  jir 

N 

(24) 


°  a -([(I  +  ^  (^)']-[(l  -  -f  *  (^)^])^ 

a  a  a  a 


=  ^icp-A:^ 

which  is  independent  of  both  <j)'  and  a. 


Approach  Velocity 

Figure  7  shows  schematically  a  plan  view  of  the  velocity  field  on  the  upstream 
surroundings  of  the  gate.  The  average  flow  velocity  through  the  e^erture  may  be 
estimated  from: 


V  = 


lull 

2aa 


nil 

2a 


1  ^ 
1  -  0.389  6; 


(25) 


where  a  is  the  coefficient  of  contraction,  and  6  =  z/h  is  the  nondimensional 
vertical  coordinate  of  the  layer  being  considered,  with  origin  at  the  gate  hinge  and 
positive  upwards,  normalized  with  respect  to  the  water  depth  at  the  gate  hinge:  h 
(i.e.,  6  e  (0,1)).  Equation  25  implies  linear  variation  of  the  coefficient  of 
contraction  from  a  =  1,  at  the  gate  hinge,  to  its  experimental  value  as  obtained  by 
Kirchhoff,  a  =  0.61  (Lamb  1945,  Prandtl  and  Tietjens  1957),  at  the  gate  ridge. 


Joukowski’s  Transformation 

Consider  the  conformal  mapping  of  circles  x^  +  in  the  complex 

Z-plane  into  the  complex  C-pIane  (C  =  ?  +  i  *  fl)  under  the  Joukowski’s 
transformation: 


+  b>0  (26) 

A  set  of  two  real  equations  is  obtained: 


,  (27a) 

1  +  - 

4i?2 
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(27b) 


y  = 


1  - 


4R^ 


If  Equations  27  are  inserted  in  the  original  circle  equation,  the  following 
family  of  ellipses  is  obtained  as  a  result  of  the  transformation: 


(R  + 


=  1 


(28) 


The  focal  distance  for  these  ellipses  is  given  by  Equation  3: 

=  (7?  +  —f  -  (R  -  —f  =  (29) 

°  47?^  ^47? 

or  ^0  =  b  (independent  of  R!). 

Therefore,  Joukowski's  transformation  maps  concentric  circles  +  y^  =  R^,  in 
the  Z-plane,  into  confocal  ellipses  in  the  C-plane,  with  eccentricities  given  by: 


^  ^  (30) 

a  ^  2R 

Note,  in  particular,  that  the  circle  [  Z  |  =  b/2  transforms  into  the  flat  focal 
segment  FjFj  (e  =1),  as  Figure  8a  shows.  On  the  other  hand,  a  circle  on  the  open 
disk  I  Z  I  <  b/2  maps  into  a  corresponding  ellipse  in  the  C-plane,  in  such  a  way 
that  the  disk  maps  into  the  whole  slit  C-plane  (i.e.,  the  C-plane  without  FjFj ).  The 
semi-infinite  open  annulus  |  Z  |  >  b/2  again  maps  into  the  whole  slit  C-plane,  as 
Figure  8b  shows. 

Velocity  Field  on  Upstream  Face  of  Gate 

Applying  the  sequence  of  complex-plane  conformal  mappings  described  in 
Figure  9,  i.e.,  inverse  Joukowski's  transform  followed  by  a  quadrant  rotation  and  a 
direct  Joukowski's  transform  (Greenberg  1978),  one  arrives  at  the  simple 
hydrodynamic  problem  involving  flow  parallel  to  a  plate,  with  obvious  solution. 
Transforming  back  to  the  original  complex  plane,  the  components  of  the  velocity 
field  are  obtained  fi'om  the  following  set  of  simultaneous  equations: 

„2  -  ,2  .  .^2.(1  , 

-  b^f  * 
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C -plane  (double) 


Rgure  8.  Joukowski’s  transformation  in  the  complex  plane  applied  to  (a)  the  circle  |z|  =  b/2,  and  (b)  a 
circle  either  internal  or  external  to  the  circle  |z|  =  b/2 


uv 


_ 

+  4x^y^ 


(31b) 


When  evaluated  on  the  vertical  projection  of  the  upstream  face  of  the  gate 
( I X  i  ^b,  y=0).  Equations  3 1  give; 


16 
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Inserting  Equation  25  in  Equation  32  gives; 


1 

u  =  ( - )  • - 

2  a  1  -  0.389  6 


This  velocity  field  is  schematically  represented  in  Figure  10.  Notice  that  u 
vanishes  for  x  =  0  ;  therefore,  the  center  line  of  the  upstream  face  of  the  gate  is 
effectively  a  stagnation  line. 


Figure  10.  Velocity  field  on  upstream  face  of  gate 
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Estimation  of  Hydrodynamic 
Pressure  Fieid  on  Gate 

The  steady-state  component  of  the  hydrodynamic  pressure  field  on  the  gate  is 
conveniently  estimated  from  Bernoulli's  energy-balance  equation  for  the 
streamline  along  the  y-axis  in  Figure  10,  from  y-oo  to  the  stagnation  point  A, 
followed  by  the  same  fonnulation  for  points  internal  to  the  segment  AB  ( |  x  |  <b). 

For  the  y-axis  streamline,  between  a  far  away  point  and  point  A,  one  has: 


Pa  = 


2 


(34) 


Using  Equation  34,  one  has  for  points  on  segment  AB; 


p(x)  =  £-(^2  -  m2) 

or,  normalizing  to  obtain  a  nondimensional  expression: 


KP) 


1 


1  -  (-) 

P 


,  IPI<1 


where  P  =  x/b . 


(35a) 


(35b) 


This  pressure  field  vanishes  for  pj  —  l/v(2,  and  becomes  suctional  toward  the 
gate  edges,  its  magnitude  increasing  without  limit  at  B.  This  happens  because  the 
magnitude  of  the  velocity  in  Equation  33  also  grows  unboundedly  as  x 
approaches  b.  On  account  of  the  viscosity  in  a  real  fluid,  the  flow  velocity  and 
the  associated  suction  will  only  attain  a  bounded  magnitude  in  the  regions  of  the 
gate  toward  the  edges.  Also,  in  an  actual  situation,  the  hydrodynamic  pressure  at 
the  very  edge  B  must  be  zero,  as  indicated  graphically  in  Figure  10.  In  an  attempt 
to  represent  these  effects  in  the  model,  the  hydrodynamic  suction  is  allowed  to 
grow  to  a  certain  value  and  then  is  assumed  to  ^proach  its  boundary  value  at  B 
(Pb  =  0)  in  a  linear  fashion  (see  Figure  11).  The  equivalent  normalized  distance 
Ap  that  defines  the  hydrodynamic  pressure  distribution  is  identified  below  by 
experimental  analysis. 


Hydrodynamic  Line  Force  on 
Vertical  Projection  of  Gate 

Integrating  the  hydrodynamic  pressure  field  over  the  gate  width,  one  obtains 
an  expression  for  the  resultant  hydrodynamic  force  on  the  vertical  projection  of 
the  gate,  per  unit  length  of  height: 
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Figure  1 1 .  Hydrodynamic  pressure  on  the  right-half  upstream  face  of  gate 


\  2  p2  -  1  2  p; .  , 

or,  normalizing  to  obtain  a  nondimensional  expression: 


fn 


<‘f,; 


=  4P,  -  In 


1  ^  Pi  2p;  -  1 

1  -  Pi  ~  1  -  P, 


’  Pi  ^  Po 


or,  in  terms  of  the  equivalent  normalized  distance  AP; 


-  Pi)  (36) 


(37) 
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(38) 


—h—  =  x(AP)  =  ^  -  ln(A  .  1) 

(Pll).b  l-(l-AP)  ^AP  ^ 

2  ^ 


Hydrodynamic  Resultant  Force 
on  Vertical  Projection  of  Gate 

Integrating  the  hydrodynamic  line  force  in  Equation  38,  the  hydrodynamic 
resultant  force  on  the  vertical  projection  of  the  gate  is  obtained  as; 


p^ 


H 


•X(AP)-| 


V^dz 


(39) 


or,  normalizing  to  obtain  a  nondimensional  expression; 

- 2  =  ^  1)'— ^•X(Ap)-sme 

pK^  2  2  (2af 

where 

Ag  =  area  of  the  upstream  face  of  the  gate 
9  =  gate  elevation  (Figure  12) 

The  location  of  the  corresponding  pressure  center  is  given  by; 


n 


or,  after  normalizing; 

6  =  (y-y-J  (-^)i-x(Ap)6rf6 

^  0 

which  leads  to; 


z  =  0.58  h 


(40) 


(41a) 


(41b) 


(42) 


which  is  a  fixed  value  for  a  given  gate  elevation  (0). 
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3  Application  of  Theory 


Reactive  Thrust  at  Support  Shaft 

Figure  12  shows  the  free-body  diagram  for  the  gate  at  elevation  0,  subjected 
to  both  static  and  hydrodynamic  loads.  W  in  the  figure  is  the  dry  weight  of  the 
gate,  whereas  W,,,  is  the  weight  of  the  water  supported  by  the  gate.  Fh  is  the 
resultant  hydrodynamic  force  on  the  vertical  projection  of  the  gate  acting  at  the 
pressure  center  defined  by  z.  P  is  the  resultant  axial  load  in  the  support  shaft 
(prop  reaction),  acting  at  an  inclination  Q.  Rh  and  Rv  are  the  horizontal  and 
vertical  components,  respectively,  of  the  resultant  reaction  at  the  support  hinge  of 
the  gate. 

The  prop  reaction,  P,  is  easily  expressed  in  terms  of  Fh  ,  by  considering  the 
rotational  equilibrium  of  the  gate; 


Once  the  prop  reaction  is  known,  the  components  of  hinge  reaction  Rh  and  Ry 
may  be  convmently  obtained  from  the  translational  equilibrium  of  the  gate; 

Rh  =  Rc  oQ  -  Ffj  (44) 

=  W  +  -  P  s  i  ffi  (45) 

and  the  function  x(AP)  can  be  identified,  if  desired,  by  inserting  Equation  40  into 
Equation  43; 


xm  =  ^si^^Q)-(^^ifjcoft 

4.7p^j^L£s  i^ 

(2af 


(46) 


The  difficulty  is  precisely  that  P  depends  on  the  hydrodynamic  resultant  Fh  , 
which  in  turn  depends  on  Ap.  This  closed  cycle  of  cause  and  effect  is  broken 
through  experimental  analysis.  Statistical  characterization  of  experimental 
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records  of  prop  reaction  (P)  in  a  small-scale  physical  model  conveniently  define 
the  hydraulic  forces  and  allow  their  prediction  for  the  prototype  gate. 


Second-Moment  Statistical  Characterization 
of  Response  Parameters 

Let  the  total  prop  reaction,  P,  consist  of  two  components:  (1)  Ps ,  the 
hydrostatic  component,  mid  (2)  Ph  ,  the  reaction  to  the  hydrodynamic  force  on  the 
gate,  Fh  (Naudascher  1991).  Thus,  from  Equation  43,  one  has  the  hydrostatic 
component: 


L  W  ^ 

P=-  cos0(-?-  +  — ) 
d  2  3 


(47) 


which  is  taken  as  deterministic  and  is  readily  computed. 

Also  from  Equation  43,  one  has  for  the  hydrodynamic  component: 


H 


(48a) 


or 


P^=  1.16. 


sin0 


sin(0  +  Q) 


H 


(48b) 


which  may  be  written  as: 


Fh  =  IKQ)-Ph 


(48c) 


where 


D/o\  •  sin(0  + 
i?(0)  ~  sm  — • — L_— 

3  sin  0 


(48d) 


Figure  13  shows  the  hydrodynamic  equilibrium  of  forces  on  the  gate.  Notice 
that,  in  general,  the  reaction  at  the  support  hinge,  R,  pulls  the  gate  with  a 
downward  vertical  component 

The  transfer  factor  in  Eqimtion  48d,  R(0),  does  not  vary  much  with  either  the 
prop  inclination  (Q)  or  the  gate  elevation  (0)  for  0  >  30°,  as  Figure  14a 
illustrates.  The  neighborhood  of  0  =  60°  is  expanded  in  Figure  14b.  Notice  that 
for  a  prop  inclination  £2  =  45°  and  gate  elevation  0  ~  60°  ,  Ph~  Fh- 
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Figure  1 3.  Hydrodynamic  equilibrium  of  forces  on  gate 


to  general,  the  expected  value,  the  standard  deviation,  and  the  coefficient  of 
variation  of  the  hydraulic  force  on  the  gate  are  respectively  given  by: 

(45 

(49 


=  v„ 


(49i 
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GATE  ELEVATION  [DEGREES] 


0.7 


in  terms  of  the  corresponding  quantities  for  the  prop  reaction. 

Previous  experience  widi  multiple  similar  records  (de  Bejar  1995;  de  B6jar 
and  Hall  1995a,b;  Fletcher  and  de  Bejar  1995)  suggests  that,  for  all  practical 
purposes,  the  distribution  of  the  peaks  of  Ph  is  approximately  normal.  Therefore, 
the  statement  can  be  made  that,  with  ~  97  percent  confidence,  the  upper  and 
lower  bounds  to  the  design  hydrodynamic  load  on  the  gate  may  be  taken  as; 


Fh  =  ±3-0^^  (50) 

applied  at  the  hydrodynamic  pressure  center: 

z  =  0.58/2  (51) 

If  the  physical  model  is  designed  according  to  Froude’s  law  with  fluid-density 
scaling  A,p  =  1,  then  the  upper  and  lower  bounds  to  the  magnitude  of  the 
hydrodynamic  force  on  the  prototype,  in  terms  of  the  same  force  in  the  model,  are 
given  by  (Novak  and  Cabelka  1981): 


•PROTOTTPE 


= 


^MODEL 


(52) 


where  Xl  is  the  scaling  factor  for  length,  and  the  hydrodynamic  forces  are  applied 
at  the  homologous  pressure  centers  (6). 

The  second-moment  characterization  of  the  flmction  x(AP)  is  readily  obtained 
from  Equation  40  (Benjamin  and  Cornell  1970); 


E[X(A^)]  = 


j^2  % 


■^sme 

<2a) 


(53a) 


o„  = 


,  7C 


!)• 


{2af 


sin  6 


Vr  =  V„  =  V  /> 
Ph  (_£) 


(53b) 


(53c) 


Given  an  experimental  record  for  Ph  in  a  physical  model,  by  assuming  the 
record  to  be  part  of  a  sample  of  an  ergodic  process,  the  moments  of  Ph  are 
estimated  from  the  average  [~  E(  Ph  )]  and  the  area  under  the  power  spectral 
density  function  [~  OpH^  ]  (Bendat  and  Piersol  1993,  Crandall  and  Mark  1963). 
Then,  the  analysis  proceeds  by  evaluating  Equations  49  through  Equations  53. 
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Alternately,  if  from  past  experience,  the  expected  value  and  the  coefficient  of 
variation  of  the  ratio  Ph/P^  may  be  assumed,  then  the  moments  of  Ph  may  be 
estimated  by: 

E(P^)=E(P^/P^)-P^  (54a) 

=  ^(Eff/P^)-E(Pjf)  (545) 

and  the  analysis  proceeds  by  evaluating  Equations  49  through  Equations  53,  as 
above. 
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4  Verification  of  Theory 


Approximate  Experimental  Verification 

The  mathematical  model  for  the  resultant  hydrodynamic  force  on  the  gate  was 
approximately  verified  from  experiments  on  a  1/15-scale  physical  model  of  a 
typical  wicket  gate  for  the  Olmsted  locks-and-dam  project.  A  model  gate  in 
closed  position  (9=  60®)  at  the  center  of  a  three-gate  aperture  faced  a  ridge-level 
upper  pool,  with  the  lower  pool  below  the  level  of  the  support  hinge  (see  the 
example  of  application  below  for  other  characteristics  of  the  physical  model). 

Under  these  hydraulic  conditions,  the  velocity  of  the  approaching  flow  at  the 
surface  was  measured  approximately  as  3.048  m/sec  (10  ^s)  at  the  center  of  the 
aperture,  from  which  Equation  25  gives  Ka=  0.072  mVsec.  Taking,  as  a  first 
approximation,  Ap~  (l-Po)/2,  Equation  38  gives  x(4P)~  0.609,  and  Equa¬ 
tion  40  provides  an  estimate  of  the  resultant  hydrodynamic  force  on  the  gate  as 
Fh  ~  27.8  N.  This  approximate  semitheoretical  estimate  should  be  compared 
with  the  average  of  Ae  directly  measured  experimental  value  of  the  same  force 
obtained  from  Equation  49a,  after  averaging  the  records  of  calibrated  strain  gages 
on  the  model  gate  shaft; 


E{F„)  =  E{P^)  =  29.15  N 


The  semitheoretical  estimate  errs  by  approximately  5  percent  on  the 
unconservative  side,  but  this  relatively  small  error  suggests  the  usefulness  of  the 
semitheoretical  model  for  engineering  applications. 


Example  of  Application 

The  procedure  for  the  statistical  characterization  of  response  parameters  is 
illustrated  with  an  example  of  application  for  the  Olmsted  locks-and-dam  project 
Figure  15  shows  a  typicd  30-second  record  for  the  total  shaft  load  on  a  1/25-scale 
physical  model  of  a  wicket  gate  at  WES,  under  ridge-level  upper  pool  and  lower 
pool  below  the  level  of  the  support  hinge.  Figure  16  shows  the  corresponding 
power  spectral  density  function.  The  hydraulic  environment  is  that  associated 
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Figure  15.  Time  history  of  shaft  load  response  (after  Chowdhury,  Hall,  and  Pesantes  (1994)) 

with  flow  through  a  three-gate  J^erture  impinging  on  a  single  central  gate  in 
closed  position  (0  =  60°).  Under  these  conditions,  the  flow  is  observed  to  be 
quasi-laminar.  Other  pertinent  quantities  for  the  physical  model  are: 

Aperture  size:  370  mm. 

Dimensions  of  the  physical  model  gate:  1 1 8mm  x  3 1 7  mm. 

Weight  of  the  physicd  model  gate  (W):  10.23  N. 

Prop  inclination  (Q):  45°. 

Kj  (identified  during  the  experiment):  0.072  mVsec. 

Unit  weight  of  water  (T  =  16  C°  ):  Yw  =  9,800  N/m^. 


Estimation  of  hydrodynamic  force  on  the  prototype  gate: 

The  weight  of  the  water  acting  directly  on  the  gate  is  given  by: 

=  \  sin(20)  =  25.15  N 

The  hydrostatic  component  of  the  prop  reaction  is  obtained  as: 

p  ^  2 -cos 60°  10.23  ^  25JS  _ 

'  sm(60  *  45)"  '2  3  ' 

which  allows  the  calculation  of  the  mean  value  of  Ph  as: 
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Figure  16.  Autospectral  density  of  shaft  load  (after  Chowdhury,  Hall,  and  Pesantes  (1994)) 

E{Pfj)  =  E{P)  -  P^=  43.15  -  14  =  29.15 

The  variance  of  the  hydrodynamic  prop  reaction,  Ph  ,  is  obtained  by 
computing  the  area  xmder  the  power  spectral  density  function  (Figure  16); 
Oph'~  0.198 


From  Figure  14b,  the  transfer  factor  is  obtained  as  R(0)  =  0.96 
Equation  49a  gives; 

E(F„)  =  R{Q)-E(Pfj)  =  0.96-29.15  =  28 
and  Equation  49b  gives; 

Op^  =  R(0)-Op^  =  0.96-0.445  =  0.43  A/' 
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The  97  percent  confidence  range  for  the  resultant  hydrodynamic  force  on  the 
physical-model  gate  is  given  by  Equation  50; 


Ph  =  Ei-Pn)  ±  3  •  =  28  ±  3  •  0.43  =  {26.1  N,  29.3  iV), 

applied  at  the  pressure  center  given  by  Equation  51; 
z  =  0.58  •  317  •  sin  60°  =  159  mm 

For  a  physical  model  designed  according  to  Froude's  law,  with  Xp  =  1,  the 
estimate  of  the  97  percent  confidence  range  for  the  resultant  hydrodynamic  force 
on  the  prototype  gate  is  given  by  Equation  52; 

^^PRcmmpE  "  ^  =  (417  kN,  458  kN) 

applied  at  z  =  4m. 
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5  Conclusions 


This  investigation  has  produced  a  mathematical  model  that  allows  an 
engineering  estimate  of  the  hydraulic  forces  on  a  prototype  wicket  gate  in 
operation  under  the  hydraulic  environment  associated  with  flow  through  a  three- 
gate  aperture  impinging  on  a  single  central  gate  in  closed  position.  This  is  one  of 
the  critical  conditions  to  be  considered  in  the  gate  design.  The  model  is  not  valid 
for  other  hydraulic  environments  involving  turbulence  in  the  back  of  the  gate. 

It  is  recommended  that  two  basic  measurements  be  taken  during  the 
experiments  with  a  small-scale  physical  model  to  completely  define  the 
theoretical  development:  (1)  the  velocity  of  the  approaching  flow  at  a  specific 
point  of  the  reservoir  nearby  the  aperture,  and  (2)  a  time-history  record  of  the 
prop  reaction  for  the  gate  in  operation  in  closed  position.  The  model  predicts 
sufficiently  close  bounds  for  the  peak  hydraulic  forces  as  to  produce  reliable 
designs. 
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